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Uncertainty in the PQCD calculation of B decays is investigated in B —> n, B —> D transition 
form factors and B — > Dtt decay amplitudes. B meson distribution amplitude dependence is studied 
by taking three kinds of distribution amplitudes so far suggested. It is found that almost same 
dependence of the form factors can be obtained irrespective of the types of the B meson distribution 
amplitudes by suitably choosing one parameter. B — > Dtt process shows the difference due to 
the distribution amplitude. The effect of the sub-leading component of the B meson distribution 
amplitude is also studied in the three processes. The numerical results of calculations with the 
sub-leading component can be well approximated by the leading order calculation with a suitable 
choice of the distribution amplitude parameters. 
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^ ■ I. INTRODUCTION 

B meson decays has been attracting much attention to check the consistency of the standard model (SM) and 
. to explore the existence of a new physics beyond the SM. Two B physics dedicated experimental facilities are con- 
structed at KEK and SLAC. The Belle and the BABAR groups have reported a lot of interesting results since their 
beginnings 1]. Many fruitful theoretical works on B physics have been made in these decades, but hadronic effects 
p\| often obscure the theoretical predictions. Li and collaborators developed the so-called PQCD method and applied it 
J> , to exclusive B meson decays as one of the approaches to tackle this issue 0. The PQCD method gives reasonable 
i predictions on B ^ Ktt^, B ttttPI and other B decays|j5|. 

I In PQCD method, a decay amplitude is obtained as a convolution of a hard part (H) and meson distribution 
1^ : amplitudes {4>k)- 

O: r 

^ ■ Amp = (l)ix H X (j)2-- - . (1) 

5 ■ 

' The hard part can in principle be perturbatively calculated in a systematic way, while the non-perturbative contribu- 
tions are incorporated into the distribution amplitudes. Major uncertainty in the PQCD calculation lies in the choice 
O il of distribution amplitudes. We need a model or a non-perturbative method like QCD based sum-rule0, Q, S to 
D , obtain the distribution amplitudes. Meson distribution amplitudes are important also in the study of B non-leptonic 
' decays with QCD factorization method^lfl] and in the calculation of the form factors with QCD based sum-rule 
^ . scheme0,l3- So far most of the PQCD calculations of B decays are given in the leading order of as and 1/AI. (A 
trial to estimate the higher order effects in as is given in ^J.) The aim of this paper is to investigate the uncertainty 
. , of the leading order PQCD calculations. Our strategy is as follows: In Sec. II, we analyze B ^ n form factors to 
' estimate the uncertainty due to the factors given below; 

1. B meson distribution amplitude: _B — > tt form factors are calculated by adopting three kinds of B meson 
distribution amplitudes proposed in the previous works [3, Il2l Il3j . The parameters of B meson distribution 
amplitudes arc fixed to accommodate with the reasonable value of the form factor at = 0. Then we vary 
these parameters to see how the value of the form factor changes. 

2. Pion distribution amplitude: We adopt the distribution amplitude given in QCD based sum-rule analysis0,0], 
and investigate the dependence on the parameters of those distribution amplitudes. (The effect of choosing 
another pion distribution amplitude is investigated in [T5|.) 

3. Hard part: The dependence on Aqcd and other renormalization group parameters are investigated. 

4. Sub-leading contributions: We estimate the 0{1/M) corrections in the hard part and the contributions from 
the sub-leading component of the B meson distribution amplitude. 
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In Sec. Ill, we analyze B D form factors. The parameters of B meson distribution amplitudes are fixed by the 
first analysis. Here, we investigate the dependence on the parameter of the D meson distribution amplitude proposed 
in jl^. In Sec. IV, we analyze B Dn decays by using the B, D and pion distribution amplitudes fixed in the 
previous analyses. The non-factorizable contribution is important in B —> Dn decays |17|. We show which B meson 
distribution amplitude gives better results by calculating the non-factorizable contribution. Sec. V is devoted to 
summary and discussions. 

II. HEAVY-TO-LIGHT FORM FACTORS 

We first analyze the i? — > tt form factors in the fast recoil region with PQCD method. We shall determine the 
parameters of B meson distribution amplitudes from the i? — > tt form factors. The i? — s- tt transition form factors 
Ff'^ and F(^'^ are defined by the matrix element. 



(Pi+P2)m- 



-9m 



where q = Pi — P2 is the lepton-pair momentum. Another equivalent definition is 

(7r(P2)|^(0)7^ii(0)|S(Pi)) = /i(g')Pi^ + /2(g')P2^ 
in which the form factors /i and /2 are related to F;^"^ and P^f^ by 



^0 



(2) 

(3) 

(4) 
(5) 



In PQCD method, the form factors pfp derived from the diagrams with one hard gluon exchange shown in 
Fig. n PQCD works best in the region with large energy transfer, i.e., with small q^. Soft contribution from the 
diagram without any hard gluon is Sudakov suppressed [ig. The formulae for the B tt form factors are given as 



q q 

FIG. 1: Leading-order contribution to F^'^ . 



fi = 167rm|C_Fr^ J dxidx2 Jbidbib2db2(l)Bixi,bi)[(j)P{x2)-(l)lix2)] 

xE{t^^^)h{xi,X2MM) , 
f2 — 16nm%CF J dxidx2 j bidbib2db2(f>B{xi,bi) 



(6) 



0^(X2)(1 -I- X277) + 2r„ - X2)<j)i{x2) - X2(I)1{X2)^ 

+2r^r^Eit'^^)h{x2,xi,b2,bi) \ , 



Eit'-^^)hixi,X2,bi,b2) 



(7) 



with 77 = 2Pi • P2/m% = 1 - the ratio r.^ = mg/ms (toq- chiral mass of pion) and the evolution factor 



Fit) = a,(Oe-^-(*)-^'(*) , 



(8) 
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where Sb and Sj^ are the Sudakov factor of kT part for B meson and pion, respectively|l85 . The hard function is 
given as 

h{xi,X2,bi,b2) = St{x2)Ko {y/xiX2rimBbi) 

X [0{bi - 62)^^0 (y/x^mBbi) Iq {y/x^m,Bb2) 

+0{b2-bl)K(^{^/x^rnBb2)h{^/x^mBbl)] , (9) 
where the factor St is the threshold resummation factor 

2i+2^r(3/2 + c) 



St{x) 



[x{l-x)r 



(10) 



V^(l + c) 

which suppresses the end-point behaviors of the meson distribution amphtudes. The hard scales t(i)'(2) g^j.g defined as 

t'^^ = max(^j?^ms, 1/fei, 1/62) , 

i'^^ = max(^xi77ms, 1/foi, 1/62) • (11) 
We investigate here the following candidates of B meson distribution amplitude: 



0«^^(a;,6) = iV«^'^a;^(l~a;)^exp 

xtub 



1 / xtub 



(l3%'^x,b) = N^'^x exp 



2 \t^KLS 
1 

1 + {bcJGN? ' 



x,b) 



iv^^^^x e{x)e{ 



2^™^ .x)Jo J:r(-2'^^^'5^ 
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(12) 
(13) 
(14) 



The first one, which we call Gaussian type, is proposed in Q. The x dependence of the second one, which we 
call exponential type, is proposed in T^l, and we take its b dependence as Lorentzian, the Fourier transform of the 
exponential function. The third one, which we call KKQT type, is obtained by solving the equations of motion under 
the approximation of neglecting 3-parton contributions 13J. Each candidate is parameterized by one parameter, loklSi 
ujgn or Akkqt- The normalization constant Nb is related to the decay constant /b through the relation 



dx 4>b{x, 0) 



(15) 



The shapes of these B meson distribution amplitude with b — Q are shown in Fig. |21 where the parameters are chosen 



KKQT 



SO that 0(0) = 0.3 as explained later. 




Gaussian 
Exponential 



0.2 0.3 0.4 0.5 



FIG. 2: 4>b{x,Q) for ujkls = 0.38, ojgjv = 0.36 and Kkkqt/Mb = 0.094, 

In Eqs. ijnji and we have included the two-parton twist-3 distribution amplitudes cf)'^ and (p^ associated with 
the pseudo-scalar and pseudo-tensor structures of the pion, respectively 6] . The contribution from the axial vector 
component (j)-^ is twist-2. The pion distribution amplitudes derived from QCD based sum rule are given as0 



07r(a;) 



if. 
h 

fn 



2J2N, 



x{l -x) 1 + a2Cll'^{\ - 2x) + 04(74/^(1 - 2x) 
1 + a2^C]l'^(\ - 2x) -I- a4pC]/^(l - 2x) 
(1 - 2x) [1 + 6a2t(10a;2 - lOx + 1)] . 



(16) 
(17) 
(18) 
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The coefficients 02 , . . . ,a2t are defined asQ 



where 



02(1 GeV) = 0.44 , 04(1 GeV) = 0.25 

a2p = 3O773 - ^pI , a4p = -(3773LJ3 + 



a2t 



pI 



5% - ]^m^3 



20 



pI 



-Pla2, 



-,2 



Ik 
Vsip) = 

t^3{p) = 



mo 

ak{p-o) , b = 11 



s{p) 



as{po) 



1 



as{p) 



as{po) 
as{p) 



81 
10 



P7Ta2) 



2(A: + l)(fc + 2) 
V3{Po) , 73^ = ^Cf + Nc, 



25 7 

^3iPo) , 73" = 'y'^f + 



_as(Aio). 

with 773(1 GeV) — 0.015, 0)3(1 GcV) = —3. The Gegenbauer polynomials are defined by 



-{3t^ - 1) , Cf\t) = i(35i4 - 30t^ + 3) , 
1(5^^-1), Cr(t) = ^(21^4 -14^2 + 1) 



(19) 
(20) 

(21) 

(22) 

(23) 
(24) 

(25) 

(26) 



(27) 




A. Numerical results 



We present the numerical results of the B ^ tt transition form factors given above. The default values for inputs 
are given as follows: 

fB = 190 MeV , U 130 MeV , Aqcd = 250 MeV {Nf = 4) , 

c = 0.3 , mo = 1.4 GeV , , . 

02 = 0.44 , a4 = 0.25 , a2p = 0.43 , ^ ^ 

a4p = 0.09 , a2t = 0.55/6 . 

The parameters in B meson distribution amplitudes are chosen as ujkls — 0.38, ojgn — 0.36 and ^.kkqt /Mb = 0.094 
so that we have F^o(O) — 0-3, which is reasonable in comparison with the sum-rule results jsj. We neglect the scale 
dependence of parameters mg and 02 , . . . , 024 in the default calculation. Its effect shall be discussed later. Monte- 
Carlo method is used to evaluate the numerical integrals. We have set the number of samples so that the statistical 
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errors in Monte-Carlo integrations may be less than 0.1%. The values of two form factors should be equal at = 0. 
The PQCD results becomes unreliable gradually at slow recoil. Our results of F^'^{q^) and F^'^{q'^) for = ~ 10 
GeV^ are shown in Table ^ and Fig. 0] It can be seen that the dependences are almost same irrespective of 
the choice of the B meson distribution amplitude. The difference is at most 4 % at = 10 GeV^. The ratio of 
each contribution from , and (p^. to the total value of (0) is given in Table ^1 It shows that the twist-3 
contribution is important as explained in p^ . 



Gaussian type with lokls = 0.38: 



(GeV^) 


0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 




0.297 0.310 0.324 0.339 0.355 0.374 0.393 0.416 0.441 0.468 0.499 
0.297 0.321 0.347 0.377 0.411 0.450 0.494 0.546 0.605 0.674 0.756 


Exponential type with o^gjv = 0.36: 


q' (GeV^) 


0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 


F^- 

pB. 


0.300 0.312 0.325 0.339 0.356 0.373 0.391 0.413 0.436 0.461 0.490 
0.300 0.323 0.349 0.378 0.412 0.449 0.492 0.542 0.599 0.665 0.743 


KKQT type with ^^^''^ - 0.094: 
niB 


q' (GeV^) 


0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 


rpBr: 


0.299 0.313 0.327 0.342 0.359 0.378 0.399 0.422 0.447 0.476 0.508 
0.299 0.324 0.351 0.381 0.415 0.456 0.501 0.554 0.615 0.686 0.770 



TABLE I: Numerical outputs of F^"" {q^) and Ff ''(g^) 



Gaussian Exponential KKQT 



>7r 


(%) 


40 


37 


40 


I 


(%) 


47 


51 


46 




(%) 


13 


12 


14 



TABLE IL The contributions from cj>^, 0? and 0*^ to the total value of F^''(0). 
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FIG. 4: The B tv form factors F_f " and F(f " as functions of (GeV^). The results by Gaussian, exponential and KKQT 
type B distribution amplitude are shown in solid, dot and dot-dashed line, respectively. 



B. Parameters in distribution amplitudes 



Each of the B meson distribution amplitudes, Eqs. (|12|) . (|13|) and (|14|l . adopted in the previous calculation has 
only one parameter ujklSj ^gn and Akkqt, respectively. The pion distribution amplitudes, Eqs. I|16f) - (I18|) . contain 
5 parameters, mj, ai, a2, ?73 and ^3. We study how the numerical outputs of the form factors at = vary with 
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these parameters. The form factor at = can be rewritten by factoring out the parameters in the pion distribution 
amphtudes as 



TUB 



+aipF''\X) + F^\X) + a2tF'^\X)\ , (30) 

where X — uikls, ^^gn or Kkkqt- The functions F^^ , F^^ . . . F^^ do not depend on the pion parameters. 

Gaussian type: In the case of Gaussian type B meson distribution amphtude, the lukls dependence can be well 
approximated within 1% precision for 0.28 < lokls 5: 0-48 by the following formulae; 

F'^^iujKLs) = 0.0623 - 0.175 {ljkls - 0.38) + 0.382 (ujkls - 0.38)^ - 0.784 {ljkls - 0.38)^ , 

F-^^iujKLs) = 0.0860 - 0.246 {lukls ~ 0.38) + 0.455 {ujkls - 0.38)^ - 0.634 {lukls " 0.38)^ , 

F^\u;kls) = 0.0784 - 0.231 {lukls - 0.38) + 0.417 (ujkls - 0.38)2 _ q 242 {lukls - 0.38f , 

FP°{ljkls) = 0.446 - 1.98 {ukls - 0.38) + 6.47 {lokls - 0.38)^ - 17.5 [ukls - 0.38)^ , 

F^'^iuKLs) = 0.153 - 0.563 [lokls - 0.38) + 1.34, [ukls - 0.38)^ - 2.05 (lukls - 0.38)^ , 

F^\ujKLs) = 0.0825 - 0.288 {lukls - 0.38) + 0.591 {ujkls ~ 0.38)^ - 0.471 {lukls - 0.38)^ , 

F^°{ujKLs) = 0.109 - 0.332 {ujkls - 0.38) + 0.635 {ujkls - 0.38)^ - 0.734 {ujkls - 0.38)^ , 

F^^iujKLs) = 0.441 - 1.37 {ujkls ~ 0.38) + 2.64 {lukls ~ 0.38)^ - 3.86 {ujkls - 0.38)^ . (31) 

The chiral mass mo — rti^/ [mu + m^) plays an important role in hadron dynamics. It gives penguin enhancement 
in B meson non-leptonic decays as pointed out in 0|. It is essential for the form factor calculation to take into 
account of the important higher- twist contributions [T^ . The chiral mass toq enters in Eqs. © and lO as 

= mo/mB and in the parameter of Gegenbauer polynomials in the pion distribution amplitudes. (See 
Eqs. H20|l - H22|) .) The dependence of the from factor is linear, while the parameter in pion distribution 
amplitudes depends linearly on l/mg. The mo dependence of a2p, a^p and a2t through can be neglected since 

= 0(10^2). The UJKLS - iTT-o dependence of F^'^{0) is shown in Fig. [51(a), where other parameters are fixed 
to the default values. The dependence on other inputs are also shown in Figs. |Sl(b)-(e). 

Exponential type: The similar calculation is done in the case of the exponential type B meson distribution ampli- 
tude. The UJGN dependence is obtained for 0.26 < ujqn < 0.46. The approximation formulae of F^^ ^ are 
given in the appendix A. The result is shown in Fig. 

KKQT type: The result in the case of the KKQT type B meson distribution amplitude is shown in Fig. |7| The 
^KKQT dependence is obtained for 0.074 < h.KKQT 1Mb < 0.114. The approximation formulae of F"*" ~ F'^'^ 
are given in the appendix A. 

These figures show that F^'^(O) depends most significantly on mo and a B meson distribution amplitude parameter 
{LOB^ ujQN or Akkqt)- The change of other parameters within reasonable range affects on F^^{0) at most 10%. 

The B meson decay constant, /b, is concerned solely with the normalization of B meson wave function in the form 
factor calculation. The normalization constant Nb enters linearly in our calculation. So if /s changes to /b + ^Ib, 
the output changes to (1 + Afs/fB) times the original value. This is also the case in the calculation of non-leptonic 
decays 

C. Intrinsic b dependence 

We investigate the uncertainty from the intrinsic b dependence of light meson distribution amplitudes, which are 
advocated by KroU et al.^9]. The b dependence of pion is taken to be the following form ,19J, 



exp 



x{l - x)b^ 
4^! 



(32) 



where a,r is the transverse size parameter of the pion. We take ~ V&t^/t: here. The variation of F^'^{0) under 
the influence of the above 6 dependence of pion distribution amplitude is shown in Table UTTl The effects of intrinsic 
b dependence is estimated to be about 10% or less. 
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0.3 0.35 0.4 0.45 



(a) 



a2 




0.3 0.35 0.4 0.45 



(b) 



m 

0.3 0.35 0.4 0.45 



(d) 



Q > 0.40 
0.35 - 0.4 

10.30 - 0.35 
0.25 - 0.30 
0.20 - 0.25 
<0.20 



a4 



m 

0.3 0.35 0.4 0.45 



(c) 



(Ok 



0)3 



(e) 



FIG. 5: Contour plots of F^^(0). (a) uikls - mo, (b) ukls - 02 (c) uikls - 04, (d) ojkls - Vs (e) ujkls - 1^3- The values of 
F^"(0) are shown by using shades as given in the sample . 




(d) (Oqn (e) (Ogn 



FIG. 6: Contour plots of F^'^(O). (a) ljgn - mo, (b) ljgn - 0-2 (c) wgat - 04, (d) ujgn - Vs (e) wgjv - W3. The values of F^'^{0) 
are shown by using shades as given in the sample . 
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a2 



I I > 0.40 

0.35-0.4 
■ 0.30-0.35 

10.25 - 0.30 
0.20-0.25 
<0.20 



0.08 0.09 0.1 0.1 



(a) 




a.4 



0,08 0.09 0,1 0.1 



(b) 



m 

0.08 0.09 0.1 0.11 

(c) ^KKC 



0.02 
0.018 
0.016 
0.014 
0.012 

0.01 



(d) 



(e) 



FIG. 7: Contour plots of ^^"(0). (a) Akkqt - mo, (b) Akkqt - 02 (c) Akkqt - 04, (d) Akkqt - V3 (e) Akkqt - i^a- 
Akkqt is given in unit of Mb- The values of F'^" {Q) are shown by using shades as given in the sample . 





Gaussian Exponential KKQT 


without b dependence 


0.297 


0.300 0.299 


0.8 X VStt/^ 


0.284 


0.285 0.286 




0.277 


0.278 0.279 


1.2 X VStt/^ 


0.269 


0.269 0.272 



TABLE III: b dependence of F^''(0) 



D. Evolution effect 



1. Gegenbauer coefficients 



The Gegenbauer coefficients in the light meson distribution amplitudes depend on the energy scale. In the PQCD 
calculation it evolves with the scale 1/6 governed by [as{l/b) / as{pLo)Y' (6 = 11 — 2iV//3), where /xq represents the 
initial scale the evolution starts with, and 7 is an anomalous dimension 0. We have investigated this evolution effect. 
Calculations are made by taking the evolution effect into account. It can be seen from Table Hvl that the effect is 
about 10%, and can be covered by the theoretical uncertainty from the variation of the Gegenbauer coefficients. 



MO (GeV) 


Gaussian Exponential KKQT 


no evolution 


0.297 


0.300 0.299 


0.5 


0.347 


0.345 0.352 


1.0 


0.294 


0.297 0.298 


1.5 


0.278 


0.282 0.280 



TABLE IV: The evolution effect on F^''(0) 
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2. Aqcd 

The QCD coupling constant as appears explicitly and implicitly through the resummation factor S in Eq. (jSJ. The 
QCD scale Aqcp determines as- Let us see how the form factor values varies depending on Aqcd- The result is 
given in TableTvl which shows that change in the form factor values is about 3% for 200 MeV < Aqcd < 300 MeV. 



Aqcd (MeV) 


Gaussian Exponential KKQT 


200 


0.299 


0.308 


0.310 


225 


0.298 


0.305 


0.306 


250 


0.297 


0.300 


0.299 


275 


0.293 


0.294 


0.294 


300 


0.289 


0.288 


0.286 



TABLE V: Aqcd dependence of F^''(0) 



3. Hard scales 



The scale of as in the expression of the form factors is determined in Eq. This choice is not unique, because 

the next leading order correction has not been calculated. There is another candidate of the scale: 

^(1) = ^(2) _ max(y^2fi^mB, y/a^niB, l/^i, 1/^2) , (33) 

The change of the value of F^^(O) under the choice of hard scale is shown in Table IVTl For reference, we also show 
the value in the case of the fixed hard scales; t^^^ = t'^' = Mb/2, Mb, 2Mb- The result shows that F^'^{0) changes 
about 10% or less depending on the choice of the form of the scale 





Gaussian Exponential KKQT 


original 






0.297 


0.300 


0.299 


Eq. 






0.288 


0.290 


0.289 


fixed t'l' 




= Ms/2 


0.286 


0.288 


0.288 


fixed t'l) 




= Mb 


0.276 


0.277 


0.277 


fixed 


= t(2) 


= 2Mb 


0.269 


0.270 


0.270 



TABLE VI: Scale choice dependence of F^"(0) 



4. Threshold resummation factors 

There is a source of theoretical uncertainty from the threshold resummation factor c in Eq. H10() . Note that this 
uncertainty, whose property differs from others like tuq, is not due to an unknown parameter, but to our parameter- 
ization. In principle, we can adopt the exact resummation result, such that no theoretical uncertainty is associated 
with it. 



E. Sub-leading contribution 

1- 0{Ahadron/mB) tcrms 

The formulae of the form factors /i and /2 are the leading order results where the terms proportional to xi 
^hadron/iTiB are ncglectcd. If we do not neglect them, the following terms are added. 

A/i = 167rm%CF J dxidx2 y &id6i52d&2</'B(a;i, 6i)a;i(?70,r - Sr^r^^) 
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xE{t^^^)hix2,xi,b2,h) , (34) 
A/2 = Wnni^CF J dxidx2 J feid&i62d&2</>B(a;i, ^-(f).,, + —^(j)^^ 

xE{t^^)h{x2,xi,b2,h) , (35) 



The numerical outputs of the above quantities are given in Table IVIII It can be fond that the sub-leading contribution 
from xi terms is about 4% of the leading value. 



Gaussian Exponential KKQT 
0.040 0.035 0.041 



TABLE VII: AF^'"(0)/i^^'"(0) 



The leading order results of PQCD calculation are obtained under the approximation of Mb — mi,. If we do not 
take this approximation the formulae become as follows; 

/i = 167rm|C_Fr^ J dxidx2 J bidbib2db24>Bixi,bi)[{l + rB)(j>!^{x2) - (l - rB)(jjlix2)] 

xE{t^^^)hB{xi,X2MM) . (36) 
/2 — l%T:m\CF j dxidx2 j bidbib2db2(j)B{xi,bi) 



Mx2){l + X2V - rs) + 2r, {{I - ^)(i - 0^2)0^2^2) - a;2<^^(x2)j 



E{t'^^^)hB{xi,X2M,h 



2) 



+2r^r.E{t<-^)h{x2,xi,b2,bi)j , (37) 

where = A/ Mb = {Mb — mi,) /Mb- The hard function Hb is given as 
hB{xi,X2.,biM) = St{x2)Ki3 {y/xiX2T]mBbi) 

X0{x2ri - 2r_B) 9{bi - 62)^^0 (^\/x>2V~^2r^mBbi^ Iq (^\/ X2'n - 2rBmBb2^ 

+9{b2 - bi)Ko (^y/x^r]~^2r^mBb2^ la (^\/x^ri~~2r^mBbi^ . (38) 

The outputs of F^'^{0) with the above formulae are given in Table IVlTll Comparing this result with the leading order 
one, we find that the effect of this approximation is about 2%. 



Gaussian Exponential KKQT 
0.302 0.305 0.304 

TABLE VIII: F^''(0) without taking Mb = nit 



2. Another component of B distribution amplitudes 

The B meson distribution amplitude in fact consists of two components |l 2j| : 

$B = + ""Bhs [MUk)+ ^-<^B(fc)] . (39) 

where v = P/mB = w+ + W- with v+ = ((u" + )/V2, 0, Ot) and u_ = (0, (u° - )/\/2, Ot). The spatial direction 
of the velocity v is taken along the third direction (w^ = v'^ = 0). By using the identity (/P + 7^3)75(1+ /)) = 0, we 
can add an arbitrary function / in the above expression; 

$B cx (/P-|-mB)75 /;_0;b] 
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= {/P + mBh5 + (1+ ^)/] 

= (/P + ms)75 [/+ + + 

- -{^ + mB)75 [(/ + 0^ + + /)+ + /)] 

In the rest frame of B meson, v'^ = 1 and other components of v vanish, so that we have 



(40) 



(41) 



where n = (0,1, Ot) and n — (1,0, Ot)- We have so far considered the contribution from the first term alone by 
choosing / = —(j)^ or / = —(j)^, and that from the rest of the terms has been neglected. Here we estimate the 
contribution from the rest of the B meson distribution amplitude components. 

A care is necessary in choosing 0^ in the rest frame of B meson where we need to distinguish "+" direction. In 
[T^ and the coordinate of the light quark in B meson is denoted as z which is on the light-cone, = z^ = 0. 



\q{z)TKmB{p)) 



-Tr 



75r— — < (t)+{v ■ z)- /z 

2 2v ■ z 



(42) 



where (j)± are the Fourier transforms of (f)± . The function 0^ is defined as the distribution amplitude associated with 
f"*", so that it becomes the leading distribution amplitude in the limit t ^ v ■ z = z~ -\-v~ z~^ — )■ oo. Since z'^z~ — 0, 
z'^ — Q and 7^ are taken in their treatment. Then the momentum of the light quark, k, should be taken along 
"+" direction, so that z • fc 7^ 0. We have taken the light quark momentum along "-" direction in the calculation of 
Eqs. ® and GJ [13. Then (f) is the leading distribution amplitude, and our choice corresponds to f — —0^. 

Let us express (j)B = f + (f)^ + 4'b^ ~ 't'ts ^ f ^^"^ 't'B = 'Pb + f- The contribution from is given in Eqs. © 
and iQ. The contribut ion from 0^ is given as 



rn 
.12 



-IGvrm^Ci? J dxidx2 J 6i(i6i&2'i&20s(a;i, &i) 



+2r^<l>PE{t^^)h{x2,xi,b2,bi) 



E{t^^y)h{xi,X2M,b2) 



(43) 



(44) 



The contributions from 0^ is given as 

/f = -16nm%CFr^ J dxidx2 J bidbib2db2(t)%{xi,bi)[(j}P{x2) - (t)Ux2)] 

x£;(t(i))/i(xi, 2:2,61,62) , 
f2 — —Idnrri^CF [ dxidx2 [ bidbib2db2(t>B{xi,bi) 

E{t^^^)h{xuX2MM) ■ 



4>tt{x2) -r^{- + X2)'f^{x2) + r^(- - 2:2)0^(2:2) 
r/ 77 



Note that the sum of contributions from (pB, 4>^b ^^^^ 0^ vanishes if 0b —4'% ~ 4'b- because 



Mx) + ^^Bix) + ^^Bix) 



= [/P + ms]75[l + ^ + -^]0B(2:) 
= mB[l+ X:]75[l+ ^\<1>b{x) = . 



(45) 



(46) 



(47) 



We need 0^ and 0^ to calculate the numerical values of these contributions. The candidates of the leading 
distribution amplitude, 0^ in this case, are already given in Eas. (|12l) - H14I) . For KKQT type distribution amplitude, 
0^ is derived in J^. (Note that the 0^ in corresponds to 0^ here.) 

.^KKQT 



6+(™)(x,6)^7Vb(2: 



rtiB 



:c)9ix)9{'^JtI£ISQL _ ^) j„ bjx{^^^^^ - x) 



UlB 



niB 



(48) 
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The X dependence of the candidate in the case of exponential type is proposed in . We add the same b dependence 
as in Eg. lfT^ . 



\ ITlB ) 



exp 



xrriB 

^GN 



1 



1 + (buGN)^ 



(49) 



As for the Gaussian type case, a candidate of (/)^ is proposed in by solving the equation of motions given in |T3| 
with (0^ + 0^)/2 — (f)g^^{x, b). Here we take (f)'g — (t)g^^{x, b), and put it into the equation of motions. The details 
are given in the Appendix B. The result is as follows; 



KLS ^KLS 



B 



m. 



exp 



1 / xniB y 

2 \ujklsJ 



+v27rm 



jErf ( 



xniB 



\V2lo 



KLS 



exp 



{m|(l 



-^i^KLsb)^ 



2^KLS} 



C 



(50) 



where the constant C is chosen so that 



iHKLS) 

■'b 



(l,fe) = 0. 



We have investigated the contributions from 0^ components under the choice of / = —<j)%- The results of F-^'^(O) 
with both (f)^ and 0^ contributions are shown in Table IIXI It can be found that ujkls = 0.45, logn = 0.42 or 
^KKQT /rriB = 0.12 gives F^^(O) = 0.3. The dependence of F^"^ with both contributions is shown in Fig|Sl The 





0.43 0.44 0.45 0.46 0.47 


F^"(0) 


0.317 0.308 0.298 0.289 0.281 






0.40 0.41 0.42 0.43 0.44 


F^"(0) 


0.321 0.310 0.300 0.291 0.282 



J^KKQT /mB 


0.10 


0.11 


0.12 


0.13 


0.14 


F^"(0) 


0.374 


0.336 


0.303 


0.274 


0.249 



TABLE IX: The value of F^'^(O) for ujkls, ujcn and K/mB by using Gaussian, exponential and KKQT type distribution 
amplitudes, respectively. 



results with the leading contribution only {ojkls — 0.38, lugn — 0.36, Akkqt/Mb — 0.094) are also shown for 
comparison. It can be seen that there is little difference for low between two kinds of calculations. Say in other 
words, the inclusion of 0^ contribution can be well approximated just by choosing a suitable value of the parameter, 
luklS: ^GN or AxKQT- We found that the difference between the two kinds of calculations is about 3% or less for 
g2 < 5 GeV^. 

For a reference we show the ratio of the contribution from the 0^ = (0^ — (/)^)/a/2 component to that from the all 
components in Table 1x1 The (f>^ component contribution is found to be about 30% or less. 





Gaussian Exponential KKQT 




0.22 0.20 0.29 



TABLE X: Ratio of the sub- leading contribution to the total one in F^'^(O) 



III. HEAVY TO HEAVY FORM FACTORS 



In this section we investigate the heavy-to-heavy form factors in the fast recoil region, concentrating on the B D 
transition. We shall determine the parameters of the D meson distribution amplitude. The B ^ D transition form 
factors are defined by the matrix elements. 



(i?(P2)|6(0)7Mc(0)|S(Fi)) = V^^^^7^^^[^+(77)(«i-f t;2)p+^-(?7)(«i-t'2)p] , 



(51) 
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FIG. 8: The B — > tt form factors F^^ as functions of (GeV^) for Gaussian type (a), exponential type (b) and KKQT type 
(c) distribution amplitudes. The results with the leading contribution only with a suitable choice of the parameter are shown 
in dashed lines, while those with both contributions are shown in solid lines. 



where rj — Pi ■ P2/ {mBTno)- The lowest-order diagrams for the B D form factors are similar to Fig. ^replacing u 
and TT by c and D, respectively. The leading-order formulae have been derived in [T6l| : 



^+ = 167rCFVrm% J dxidx2 j 6id6i62'i&20s(a;i, 6i)0_d(x2) 

£;^(t(i))/i^(a;i, a;2, 61, 62) + ri;^(t('')/i^(x2, xi, 62, 61) 



X 

e- = 



where the color factor Cp =4/3 and r = m£i/mB- The functions E (t) and /i (xi, X2, 61, 62) are defined as 

= as{t)exp[~SB{t)-SD{t)], 



(52) 
(53) 

(54) 



/i^(a;i, 2:2,61,62) = Kq ( XiX2rri+mBbi ) St{x2) 



0{bi - b2)Ko yy/x^rq+TriBbij Iq y^/x^^rij+rnBbs 
+9{b2 - 61)1^0 (^\/x2rri+mBb2^ Iq (^\/x^rr]+ tob6i^ 



(55) 



where ri~^ — rj + yjrf — \. The definitions of the hard scales t'-"^'^) are as follows, 



^(1) = max(\/x2r77+m_B, 1/61, 1/62) , 
^(2) — max(\/xir77+mB, 1/61, 1/62) . 

For numerical estimation, we use the model of D meson distribution amplitude adopted in |16|| . 

4>d{x) = Ndx{\ - a;)[l + CoiX - 2x)] , 



(56) 
(57) 
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where Cd is the D meson distribution amphtude parameter. The normahzation constant Nd is found to be Sfn/V^^ 
by using the relation 



/ 



dx(j)£){x) 



(58) 



A. Numerical Results 



Here we investigate the distribution amphtude dependence of the B ~f D form factor. The inputs for B meson are 
same as the case oi B ^ tt form factor. The D meson distribution amphtude has only one parameter Cd- We take 
/d = 240 MeV, and other parameters are same as the B it case except for the threshold resummation parameter c 
which is taken to be 0.35 in B ^ D transition flGj. The D meson distribution amplitude H57|l is decomposed into two 
parts; 

<j)Dix) = q^"Dix)+CD^Dix) (59) 

where (f'^ix) = Ndx{1 — x) and ^^^(a;) = Ndx{1 — x){l — 2x). The contributions from (j)'^ and at ?7 = 1.58 (near 
maximal recoil) are shown in Table IXll for 3 types of the B meson distribution amplitudes. It can be seen that the 
value of ^+ varies about 4% under the 10 % change of Cd- We fix Cd to be 1.5 so that the value of ^+ agrees with 
the experimental data, ^+(1.58) ~ 0.6. The rj dependence of ^+ is shown in Fig. 9 for Cd = 1-5. The results shows 
that the B meson distribution amplitude dependence of f+ is less than 5%. 



contribution 


Gaussian Exponential KKQT 


total {Cd = 1.5) 


0.331 0.360 0.334 
0.167 0.167 0.170 
0.582 0.610 0.589 



TABLE XI: Contribution to ^+ from and 0}; 




0.4- 
0.2 

1.35 1.4 1.45 1.5 1.55 1.6 x\ 



FIG. 9: The B D form factor 5+ as a function of — vb ■ vd- The results by Gaussian, exponential and KKQT type B 
distribution amplitude are shown in solid, dot and dot-dashed line, respectively. 



B. Another component of B distribution amplitudes 

Following Sec. Ill E 21 we investigate the contributions from the another component of the B distribution amplitude 
in the B ^ D form factor. The contribution from 0^ and </>^ are given as 

= -16TTCFVrm% J dxidx2 J bidhih2dh24>'%{xi,hi)(j)D{x2,b2) 

xE^{t^^'))rh{x2,xi,b2,bi) , (60) 
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a = 



-16T:CFVrm% J dxidx2 J bidbib2db2<j)B (a^i, &i)0r>(a;2, ^2) 



and 4>^ vanishes if 0s 



as in the case of B 



The sum of contributions from (j)B, 
investigated the contributions from 

^kkqt/ttlb = 0.12, which is obtained in B ^ tt analysis. The resuhs of ^+(1.58) with both 



(61) 
(62) 

We have 



components under the choice of / = —(jjg and ujkls ~ 0.45, uogn — 0.42, 



6^ and 



contributions 

are shown in Table IXlIl It can be found that Cd — 0.6 gives ^+(1.58) ~ 0.6. The difference due to the choice of 
the B distribution amplitude becomes about 16% or less here. The (jfj^ component contribution is not numerically 
sub-leading in the case of KKQT type distribution amplitude. 

The ^+ value at -q = 1.58 changes 5~8% by the inclusion of 0^ contribution as seen by comparing the results 
given in Tables IXll and IXIll If a suitable value oi Go is taken for each B distribution amplitudes, we can reduce the 
difference. The suitable choice is Co = 0.74, 0.77 and 0.40 for Gaussian, exponential and KKQT, respectively. The 
•q dependence of f+ with both contributions with the suitable value of Cd is shown in Fig^| The results with the 
leading contribution only {lurls ~ 0.38, uign — 0.36, Akkqt/Mb — 0.094, Cd — 1.5) are also shown for comparison. 
It can be seen that there is little difference for 1.3 < ij < 1.58 between two kinds of calculations. The inclusion of the 
sub-leading contribution can be well approximated just by choosing a suitable value of the parameter, Cd as in the 
case of B vr. We found that the difference between the two kinds of calculations is about 2% or less. 



contribution 


Gaussian Exponential KKQT 


C+(0"i>) 


0.246 


0.277 


0.220 




0.176 


0.165 


0.265 




0.124 


0.130 


0.111 




0.093 


0.089 


0.153 


total {Cd = 0.6) 


0.552 


0.573 


0.643 



TABLE XII: Contribution to ^+ and from (f)%, and 



IV. B ^ Dn 

The decay rates of B — > Dir is given as 

= ^Gm,\'\V^,\'^\M^ , (63) 

where r = mD/mB- The indices, i = 1, 2, and 3, denote the modes B^ — > D+tt", B^ — > D^tt^ and B~ — > D'^tt~ 
respectively. The decay amplitudes Mi are written as 17j 

Ml = A^xt + /i3<^cxc + Mc^t + Mc^c , (64) 
M2 = l=[fD^int - ./sCcxc + Mint - A^cxc] (65) 

V2 

M3 = McKt + /D^int + Mc^t + Mint ■ (66) 

The factor ^oxt denotes the factorizable external M^-emission contributions. The factors ^int and ^cxc represent the 
factorizable internal H^-emission and Il^-exchange contributions, respectively. The amplitudes A^oxt, A^int, and A^cxc 
are the non-factorizable external PF-emission, internal M^-emission, and VF-exchange contributions, respectively. The 
factor ^cxt (6nt) is obtained by the convolution between the Wilson coefficients and B ^ D (tt) form factor. The 
leading formulae of these expressions are given in [T^ . They are summarized with the (j)^ and 0^ contributions in 
the Appendix C. 

Let us first show the leading order calculation for each B distribution amplitude without n and n contributions. The 
parameters are the same in the cases of the form factor calculations, {ujkls = 0.38, lugn = 0.36, Akkqt/Mb = 0.094 
and Cd = 1.5) The result is shown in Table Ixllll Our result of Gaussian case is slightly different from that given in 
p7|. It is partly due to the choice of the parameters and partly due to the change of anomalous dimension adopted 
in the Sudakov factor [23. We should look at the ratios between branching ratios rather than the magnitudes of the 
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. 4 
0.2 



1.35 1^4 1.45 1.5 1.55 1.6 T| 

(a) 




0.4 
0.2 

1.35 1.4 1.45 1.5 1.55 1.6 T] 

(0 

FIG. 10: The B ^ D form factors ^+ as functions of t) for Gaussian type {Co ~ 0.74) (a), exponential type {Cd = 0.77) (b) 
and KKQT type {Cd = 0.40) (c) distribution amplitudes. The results with the leading contribution only with a suitable choice 
of the parameter are shown in dashed lines, while those with both contributions are shown in solid lines. 

branching ratios since there is uncertainty in the decay constants of heavy mesons which gives overall normahzation 
of the distribution amphtudes. BR(£)+7r^) is slightly larger than the experimental data, while BR(Z?°7r°) is slightly 
smaller than that. 



decay mode 


Gaussian Exponential KKQT Exp. 


D°7v- 

D+TT- 
U TT 


5.3(1.0) 5.2(1.0) 5.2(1.0) 4.98 ± 0.29 (1.0) 

3.2 (0.60) 3.7 (0.71) 3.3 (0.63) 2.76 ± 0.25 (0.55 ± 0.06) 

0.18 (0.034) 0.11 (0.021) 0.20 (0.039) 0.291 ± 0.028 (0.058 ± 0.007) 



TABLE XIIL The branching ratios of B ^ Dtt decay modes in the unit of 10 . The number in the parenthesis is the ratio 
to BR(D''7r~). The experimental data is from |23 |. 

There is a cancellation between the Wilson coefficients Ci and C2 in the evaluation of ai{t) = Ci{t) + C2{t)/Nc 
which enters in ^int and ^cxc as can be seen in Fig llll ai{t) almost vanishes around t ~ Mb/'2'- The contribution 
from ^int is numerically significant in the B'^ — )■ D^t:^ decay amplitude. (The contribution from ^oxc is negligible.) 
For reference, we show how the branching ratios change if we adopt the fixed scale for the evaluation of the Wilson 
coefficients and as in Table IXTVl The result shows that the choice of the scale t can give large uncertainty in B — > Dtt. 

Let us estimate the 0^ and 0^ contributions. The formulae of B ^ Dtt amplitudes in PQCD calculation are 
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1.5 \ 

1 ■ 
0.5- 

-0.5 / 
-1 / 



■^4 6 



FIG. 11: Scale dependence of the Wilson coefficients. 



decay mode 


Gaussian Exponential KKQT 


fixed t — 


Mb/2 


D°n- 
D+TT- 


6.8 (1.0) 6.5 (1.0) 6.7 (1.0) 
2.6. (0.38) 2.8 (0.43) 2.6(0.39) 
0.44 (0.065) 0.34 (0.052) 0.46 (0.069) 


fixed t — 


Mb 


D+TT- 
D%« 


7.4 (1.0) 7.1 (1.0) 7.3(1.0) 
2.2 (0.30) 2.4 (0.34) 2.2 (0.30) 
0.66 (0.089) 0.53 (0.075) 0.69 (0.095) 


fixed t — 


2Mb 


D^n- 
D+n- 

L) TV 


8.0 (1.0) 7.6 (1.0) 7.9 (1.0) 
2.0 (0.25) 2.2 (0.29) 2.0 (0.25) 
0.87 (0.11) 0.72 (0.95) 0.89 (0.11) 



TABLE XIV: The branching ratios of B 
parenthesis is the ratio to BR(D''7r^). 



Dn decay modes for fixed RGE scale in the unit of 10 ^. The number in the 



obtained under the following choice of the light quark momenta in B meson [Tt^: 

xiMb 



ki 



V2 

~7r 



(1,0,0t) + A:iT for 6nt, Mint, 
(0,l,OT) + fciT for others. 



(67) 
(68) 



Then we should take the leading B meson distribution amplitude as 0^ in ^int and Mint, while (/)^ in others. (Remind 
the discussion given in Sec. HIE 21 ) The parameters of the distribution amplitudes are taken as ujkls = 0.45, 
i^GN — 0.42, Akkqt/Mb — 0.12 and Cd = 0.6. The ratios of the cj)^ and 0^ contribution to the total one in 



each component of the decay amplitude are shown in Table IXVI KeM^xc receives large contributions from cj}^ and 
But its magnitude is far smaller than those of ^ext, Cint and Mint, so that the effect is not significant in the 
total amplitudes Eas. H64|l - H66() . The <j>g and ipg contribution to Mext vanishes since X3 and (1 — X3) terms cancels 
in M2xf The branching ratios calculated with this set of parameters are given in Table IXVH BR(_D+7r~) gets lower 
and approaches the experimental value from the view point of the ratio. BR(D°7r'') becomes larger except for KKQT 
type case, which is a good tendency to realize the experimental values. The Gaussian type distribution amplitude 
becomes the best candidate here. 





Gaussian Exponential KKQT 




0.42 


0.37 


0.55 




0.16 


0.13 


0.23 


ReMe^t 


0.0 


0.0 


0.0 


ImMext 


0.0 


0.0 


0.0 


ReMi„t 


0.37 


0.38 


0.50 


ImMint 


-0.08 


-0.04 


-0.05 




-1.13 


-0.85 


-1.25 


ImMexc 


0.16 


0.20 


0.25 



TABLE XV: The ratios, {(p^ and (j!>J contribution) /(total), in each amplitude with Cd = 0.6. 
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decay mode 


Gaussian Exponential KKQT 


D+Tv- 

U TV 


6.5 (1.0) 6.2 (1.0) 8.4 (1.0) 
3.2 (0.50) 3.5 (0.57) 4.5 (0.54) 
0.25 (0.038) 0.17 (0.027) 0.27 (0.032) 



TABLE XVI: The branching ratios of each decay modes in the unit of 10 with Co = 0.6. The number in the parenthesis is 
the ratio to BR{D°tv~). 



Next we take the parameters as Cd = 0.74, 0.77 and 0.4 for Gaussian, exponential and KKQT type distribution 
amplitudes, respectively as done in the case of the B D form factor calculation. The ratios of the decay amplitude 
with (f>]j and (f/g contributions to that of the leading calculation adopted in obtaining Table IXIIII are shown in 
Table IXVIII It can be found that the leading calculation gives a good approximation with the uncertainty about 
20%. The branching ratios in this calculation are given in Table IXVlIII This result also shows a good tendency to 
approach the experimental value in comparison with the result of the leading calculation given in Table IXIIII The 
KKQT type distribution amplitude becomes the best candidate in this case. 



Gaussian Exponential KKQT 



A(L»V-) 
A{D+-K-) 
4(75%°) 



1.13 (4.3°) 1.13 (2.8°) 1.20 (5.5°) 
1.06 (-0.6°) 1.04 (-1.1°) 1.08 (-0.5°) 
1.13 (22°) 1.18 (27°) 1.22 (29°) 



TABLE XVII: The ratios of the decay amplitudes, (calculation with n and n contribution)/ (leading calculation given in 
Table XIII) in each decay mode with Cd = 0.74, 0.77 and 0.4 for Gaussian, exponential and KKQT type distribution amplitudes, 
respectively. The number in the parenthesis is the phase. 



decay mode 


Gaussian Exponential KKQT 


D°n- 
D+Tv- 

U TV 


6.9 (1.0) 6.7 (1.0) 7.7 (1.0) 
3.6 (0.52) 4.0 (0.60) 3.8 (0.50) 
0.23 (0.034) 0.15 (0.022) 0.30 (0.040) 



TABLE XVIII: The branching ratios of each decay modes in the unit of 10 ^ with Cd = 0.74, 0.77 and 0.4 for Gaussian, 
exponential and KKQT type distribution amplitudes, respectively. The number in the parenthesis is the ratio to BR(_D''7r~). 



V. SUMMARY AND DISCUSSIONS 



We have analyzed the uncertainty in the PQCD calculations of B ^ tt, B ^ D form factors and B Dn decay 
rates. The sources of uncertainty in B ^ n (D) form factors are summarized in Table reftbl-errs. The uncertainty 
in the perturbative hard part is less than 10%. The major source of uncertainty comes from the meson distribution 
amplitudes. The meson distribution amplitude is a non-perturbative quantity, so that we need a model or a non- 
perturbative method to evaluate it. The leading PQCD results varies 10~ 30 % by changing the parameters in the 
meson distribution amplitudes. The uncertainty from the RGE scale choice is small in the form factor, while it is 
large due to subtle cancellation between Wilson coefhcients in i3 ^ Dtt. 

Here we have tried three kinds of the B meson distribution amplitudes. Two of them are models and one is derived 
from the equations of motion under the neglection of 3-parton contributions. It is surprising that the three types of 
B meson distribution amplitudes give almost same PQCD results of i? — > tt (D) form factors by suitably choosing 
their parameters although the functional forms of them are rather different with one another. The non-factorizable 
contributions in non-leptonic B decays can be of help to discriminate the B meson distribution amplitudes. 

The formally sub-leading component of the B meson distribution amplitude gives significant contributions to B 
decays. This component is neglected in many of the previous PQCD calculations. But the leading B meson distribution 
amplitude alone can give a good approximation if we suitably choose the parameters. The difference can be reduced 
to be a few % for the form factors, and about 20 % for i? ^ Dtt amplitudes with a suitable parameter choice. So the 
results of the previous PQCD studies are still useful. 
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mode 


source 


uncertainty 




wb, mo 


30% 




02 , . . . , a2t 


10% 




b dependence of pion 


10% 






normalization 




evolution efltects 


10% 




Aqcd 


3% 




choice of the hard scale 


10% 




xi terms 


4% 




another B distribution amplitude 


20~30 % *^ 


B-^ D 


Cd 


4% 




another B distribution amplitude 


40~60 % 



TABLE XIX: The uncertainty from each source. *' As for the uncertainty from the another B distribution amplitude, the 
effects of the inclusion of another B distribution amplitude can be well approximated by a single B distribution amplitude with 
a suitable value of the parameter as shown in Figs. 8 and 10. 
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APPENDIX A: APPROXIMATION FORMULAE OF B ^ n FORM FACTORS 

In the case of exponential type B meson distribution amplitude the ujgn dependence can be well approximated by 
the following formulae for 0.26 < cogn < 0.46; 

F^°(wGw) = 0.0597 - 0.184 (wgat - 0.36) + 0.459 {ujgn - 0.36)^ - 1.04 {logn - 0.36)^ , 

F'^'^{ujgn) = 0.0780 - 0.235 {logn - 0.36) + 0.537 {lugn - 0.36)^ - 1.06 {logn - OMf , 

F'^^iuiGN) = 0.0702 - 0.212 {lugn - 0.36) + 0.485 {ujgn - 0.36f - 0.701 {lugn - OMf , 

F^°{lugn) = 0.507 - 2.38 {lugn - 0.36) + 8.60 {lugn - OMf - 24.7 {lugn - 0.36)^ , 

F^'^{lugn) = 0.144 - 0.537 {lugn - 0.36) + 1.52, {lugn - 0.36)^ - 3.55 {lugn - 0.36)^ , 

F^^{lugn) = 0.0748 - 0.259 {lugn - 0.36) + 0.649 {ujgn - 0.36)^ - 0.994 {lugn - OMf , 

F'^°{lvgn) = 0.0986 - 0.304 {ljgn - 0.36) + 0.711 {ujgn - 0.36)^ - 1.37 {lvgn - 0-36f , 

F'^^{lugn) = 0.394 - 1.23 {lugn - 0.36) + 2.84 {lugn - 0.36)^ - 5.09 {lugn - 0.36)^^ . (Al) 

In the case of KKQT type B meson distribution amplitude the Akkqt dependence can be well approximated by 
the following formulae for 0.074 < Akkqt /Mb < 0.114; 

F"^" {Akkqt) = 0.0604 - 0.662 {Akkqt/Mb - 0.094) + 5.55 {Akkqt/ Mb - 0.094)2 _ ^g g {Akkqt/ Mb - 0.094)3 
F^^{Akkqt) = 0.0876 - 0.946 {Akkqt/Mb - 0.094) + 6.44 {Akkqt/Mb - 0.094)^ - 56.2 {Akkqt/Mb - 0.094)^ 
F^^{Akkqt) = 0.0817 - 0.887 {Akkqt/Mb - 0.094) + 5.61 {Akkqt/Mb - 0.094)^ - 27.3 {Akkqt/Mb - 0.094)^ 
F'''' {Akkqt) = 0.450 - 7.45 {Akkqt/Mb - 0.094) + 90.9 {Akkqt/Mb - 0.094)2 _ ^^^gg {K^j^^t/Mb - 0.094)3 ^ 
F'"'{Akkqt) = 0.157 - 2.12 {Akkqt/Mb - 0.094) + 19.3, {Akkqt/Mb - 0.094)2 _ 293 {Akkqt/Mb - 0.094)^ , 
F^\Akkqt) = 0.0864 - 1.07 {Akkqt/Mb - 0.094) + 7.57 {Akkqt/Mb - 0.094)2 _ ^gQ ^Akkqt/Mb - 0.094)^ , 
F^'\Akkqt) = 0.113 - 1.27 {Akkqt/Mb ~ 0.094) + 8.19 {Akkqt/Mb - 0.094)2 _ g^ 7 (A^^^^/Mb - 0.094)^ , 
F^^{Akkqt) = 0.465 - 5.35 {Akkqt/Mb - 0.094) + 33.5 {Akkqt/Mb - 0.094)2 _ 300 {Akkqt/Mb - 0.094)^ . 
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APPENDIX B: 0" IN GAUSSIAN TYPE DISTRIBUTION AMPLITUDE 



The equations of motion for cj)'^ and cpg are given with the approximation of neglecting 3-parton contributions 



as Hi 



+x(f)g (x) = , 



2A 



X ] (j)%{x) + x(j)g{x) = 0, 

niB 



(Bl) 
(B2) 



where A — niB — mi, is the hadronic scale of HQET. By solving Ea. (jBl|) with (p^ = 't's^^ obtain Ea. (|50|l . As for 



Ea. (jB2|) the left hand side dons not necessary vanishes, but its value is less than 10% of (t>]j{0) for A/ms ~ 0.1. 



APPENDIX C: Dn FORMUALE 



The contributions to ^ext are given as 

^ext 



.1 nl/A 

IQ-KCF^/rm^ I dxidx2 / 6id6i62d620s(a;i, 6i)(/'z)(a;2, ^2) 
Jo Jo 

xasit)a2{t) expi-Seit) - Soit)] 

X [h{xi,X2,bi,b2) +rh{x2,xi,b2,bi)] , 
.1 .i/A 

167rCFVrm% / dxidx2 / bidbib2db2(j)'%{xi,bi)cj)D[x2.h2) 
Jo Jo 

xas{t)a2{t) exp[~SB{t) - Soit)] a;2, 61, 62)] , 

/•I /-i/A 

16TrCF\/rm% / dxidx2 / bidbib2db24)%{xi,bi)(t>D{x2,b2) 
Jo Jo 

xas[t)a2{t) eM-Ssit) - Suit)] 
X [~rh{x2,xi,b2,bi)\ , 



where 02 = C2 + Ci/Nc, and Ci,2 are the Wilson coefRcients. 
The contributions to ^int are given as, 



.1 .l/A 

167rC_F -\/Fm^ / dxidxs / bidbib^dbscjjBixijbi) 
Jq Jo 

Xas(tint)ai(iint) exp[-5B(iint) - S-^{tint)] 

[{1 + X3)^^{X3) + rail - 2a;3)(0^(x3) + (/)Ua;3))] /i(a;i, a;3(l - r^), 61, 63) 



+2ro0P(x3)/i(a;3, xi(l - r^), 63, 5i) 



.1 .l/A 

IGttCf V^^TT-B / dxidxs / &i(i&i63(i630^(x'i, 61) 
Jo Jo 

Xas(iint)ai(iint) exp[~5B(tint) - 'S',r(iint)] 

x[-a:;30x(a;3) - 7-o(l - 2^3)('/'^(a;3) + '?!'5r(a;3))]/i(a;i, 2:3(1 - r^),6i,63) , 



-2ro(/>^(a;3)/i(a;3,xi(l - r^), 63,61) 



.1 .l/A 

167rCi?\/rm^ / dxidxs / bidbib^db^cjigixi^bi) 
Jo Jo 

Xas(iint)ai(iint) exp[-5_B(ti„t) - 5',r(tint)] 

x[-0^(a;3) +roa;3(0^(x3) + (/)5r(x3))]ft(xi, 2:3(1 - r2),6i,63) 



(CI) 



(C2) 



(C3) 
(C4) 



(C5) 
(C6) 



(C7) 



(C8) 



where ai — Ci + C2/Nc- 



The form factor ^exc is written as 

/•I rl/A 
Cexc = 167rCi^V^m| / dx2dX3 / &2<i&2^3rf&3'/'D(a;2, 62) 

Jo Jo 

Xas(iexc)ai(iexc) exp[-S'_D(texc) " •S'^lW)] 

-a;3(?i^(a;3)/ia(a;2,a;3(l - r^),62,&3) 

The -B distribution ampHtude does not enter in ^excj so there is no (j)^^^^ contribution here. 
For the non-factorizable ampHtudes, their expressions are 



32TTV2NCFVrm% / [da;] / bidbib3db3(j)B{xi,bi)(j)jj(,){x2,bi)(j)^{x3) 
Jo Jo 



C (t ) 

xas(4)— ^ exp[-S'(ib)| 

61=62] 



xshl '{xi,bi) - (1 - a;3 + X2)hl '{xi,bi) 



327rV2iVC_FVrm| / [dx] / 6i(ibi63(i63(/)5(a;i, (x2, 6i)(/)7r(a:;3) 
Jo Jo 

X(^s{tb)^^ exp[-5(i6)|6,=6,] 



-X3/i,^^^(a;i,&i) + (1 - X3)hl^'{xi,bi) 

= 327rV2NCFVrml / ^a;] / bidbib3db3(t>B{xi,h 
Jo Jo 

C (t ) 



)4>D('){x2,bi)(j)Tr{x3) 



X2hf\xi,bi) 



Mr 



/■I /-l/A 

= Z2'K\/2NCF\/rm\ / [da;] / b\dbib2db24>B{xi,bi)(j)jjM{x2,b2)(l)-K{x3) 
Jo Jo 



Ml, 



X a« (i d ) ^^^^ exp [- 5 d) 1 63 =61 ] 



(-a;2 - X3)h^];\xi,bi) + (1 - a;2)/i^^^(a;i, 6^) 
/■I /-l/A 

2>2'K\rlNCF\frmB \ [da;] / b\db\b2db2(^^{x\,b\)(^F>M{x2,b2)<^T,{x3) 
Jo Jo 

C2{t6) 



xasitd)- 



N 



lo Jo 

■exp[-5(id)|63=6i] 



X3h';l\xi,bi) 



= 32'7rV2NCFVrm] 



: f\dx] f 

Jo Jo 



1/A 



bidbib2db2(j)B{xi,bl)(l)DM{X2,b2)(pn{X3) 



X as {td ) ^^^^ exp [- 5'(i d) 1 63 =61 ] 
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X2h^a\xi,bi) - (1 - X2)h^^' {xi,bi) 



(2)/ 



(C15) 



Ml 



Ml 



32nV2NCFV^'m% / [dx] / bidbih2dh2(t)B{xi,bi)(j)D{x2,b2)4>T,{xs) 
Jo Jo 

C it ] 



xsh'-pixi, bi) - X2hy'{xi,bi) 



(2), 



32nV2NCF\/rm% / I'^^;] / bidbib2db2(t)Bixi,bi)4>D{x2,b2)(t)n{x3) 
Jo Jo 

C (t ) 

xoisitf) eM-Sitf)\b,=b,] 
-x^h'-p {xi.bi) 

327rV27VCi.Vfm| / [dx] / bidbib2db2(t>B{xiM)4>D{x2M)(l>Ax3) 
Jo Jo 

X2h^f\xi, bi) 



(C16) 



(C17) 



(C18) 



The definitions of the functions, ha, hb and so on are given in Note that the sum of contributions from (/i^ , (j)'^ 



and 0^ vanishes if (/)_b = ^^"^ = 0^ 



[1] [http : //belle . kek . jp/| 

http : //www-public . slac . Stanford. edu/babar/ 
[2] H-n. Li and G. Sterman, Nucl. Phys. B381, 129 (1992); H-n. Li and H.L. Yu, Phys. Rev. Lett. 74, 4388 (1995); Phys. 

Lett. B 353, 301 (1995); Phys. Rev. D 53, 2480 (1996); C.H. Chang and H-n. Li, Phys. Rev. D 55, 5577 (1997); T.W. Yeh 

and H-n. Li, Phys. Rev. D 56, 1615 (1997); H.Y. Cheng, H-n. Li, and K.C. Yang, Phys. Rev. D 60, 094005 (1999). 
[3] Y.Y. Keum, H-n. Li and A.L Sanda, Phys. Rev. D 63, 054008 (2001); Y.Y. Keum and H-n. Li, Phys. Rev. D 63, 074006 

(2001). 

[4] C. D. Lh, K. Ukai, and M. Z. Yang, Phys. Rev. D 63, 074009 (2001). 

[5] See, for example, H-n. Li, in Proceedings of the 32nd International Conference on High Energy Physics (ICHEP 04), 
Beijing, China, 2004, edited by H. Chen, D. Du, W. Li and C. Lu (World Scientific, Singapore, 2004), vol. 2, p. 1101 
(erpint hep-ph/0408232 ) and references therein. 

[6] P. Ball, V.M. Braun, Y. Koike, and K. Tanaka, Nucl. Phys. B529, 323 (1998); P. Ball and V. M. Braun, 'hep-ph/98 08229l 

[7] V.M. Braun and LE. Filyanov, Z. Phys. C 48, 239 (1990); P. BaU, J. High Energy Phys. 01, 010 (1999). 

[8] P. BaU and R. Zwicky, Phys.Rev. D 71, 014015 (2005). 

[9] P. BaU and M. Boglione, Phys. Rev. D 68, 094006 (2003). 
[10] M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); Nucl. Phys. B591, 313 (2000). 
[11] S. Mishima and A.L Sanda, Prog.Theor.Phys. 110, 549-561 (2003). 
[12] A.G. Grozin and M. Neubert, Phys. Rev. D 55, 272 (1997). 

[13] H. Kawamura, J. Kodaira, C.F. Qiao, and K. Tanaka, Phys. Lett. B 523, 111 (2001); Erratum-ibid. 536, 344 (2002); Mod. 
Phys. Lett. A 18, 799 (2003). 

[14] Recently P. Ball, V.M. Braun and A. Lentz gave an update of pion distrubution amplitude parameters in hep-ph/0603063] 
(2006). They quoted, for example, a2 = 0.25 ± 0.15. The change of the central value of these parameters gives about 10% 
effect on the form factor value as shown in Sec. HB. 

[15] T. Huang and X-G. Wu, Phys. Rev. D 71, 034018 (2005). 

[16] T. Kurimoto, H-n. Li, and A.L Sanda, Phys. Rev. D 67, 054028 (2003). 

[17] Y.Y. Keum, T. Kurimoto, H-n. Li, CD. Lu, and A.L Sanda, Phys. Rev. D 69, 094018 (2004). 
[18] T. Kurimoto, H-n. Li, and A.L Sanda, Phys. Rev. D 65, 014007 (2002). 



23 



[19] R. Jakob and P. KroU, Phys. Lett. B 315, 463 (1993); P. KroU and M. Raulfs, Phys. Lett. B 387, 848 (1996); M. Diehl, 

P. KroU and C. Vogt, eConf C010430;T21 ( hep-ph/0105342 1 (2001). 
[20] Z.T. Wei and M.Z. Yang, Nucl. Phys. B642, 263 (2002). 
[21] H-n. Li and H-S. Liao, Phys. Rev. D 70, 074030 (2004). 

[22] Particle Data Group, S. Eidelman et al., Phys. Lett. B 592, 1 (2004) and 2005 partial update for edition 2006. 



